Importance of Correlation Effects on Magnetic Anisotropy in Fe and Ni 
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We calculate magnetic anisotropy energy of Fe and Ni by taking into account the effects of strong 
electronic correlations, spin-orbit coupling, and non-collinearity of intra-atomic magnetization. The 
LDA+U method is used and its equivalence to dynamical mean-field theory in the static limit is 
emphasized. Both experimental magnitude of MAE and direction of magnetization are predicted 
correctly near (7 = 4 eV for Ni and U = 3.5 eV for Fe. Correlations modify one-electron spectra 
which are now in better agreement with experiments. 



PACS numbers: 71.15-Mb 71.15.Rf 71.27.+a 75.30. Gw 75.40.Mg 



The calculation of the magneto-crystalline anisotropy 
energy (MAE) of magnetic materials containing 

transition-metal elements from first principles electronic 
structure calculations is a long-standing problem. The 
MAE is defined as the difference of total energies with 
the orientations of magnetization pointing in different, 
e.g., (001) and (111), crystalline axis. The difference is 
not zero because of spin-orbit effect, which couples the 
magnetization to the lattice, and determines the direc- 
tion of magnetization, called the easy axis. 

Being a ground state property, the MAE should 
be accessible in principle via density functional theory 
(DFT) Despite the primary difficulty related to the 

smallncss of MAE (~ 1 /xeV/atom), great efforts to com- 
pute the quantity with advanced total energy methods 
based on local density approximation (LDA) combined 
with the development of faster computers, have seen suc- 
cess in predicting its correct orders of magnitudes [||qT2) . 
However, the correct easy axis of Ni has not been pre- 
dicted by this method and the fundamental problem of 
understanding MAE is still open. 

A great amount of work has been done to understand 
what is the difficulty in predicting the correct axis for Ni. 
Parameters within the LDA calculation have been varied 
to capture physical effects which might not be correctly 
described. These include (i) scaling spin-orbit coupling 
in order to enlarge its effect on the MAE p|JlO||, (ii) cal- 
culating torque to avoid comparing large numbers of en- 
ergy poll , (iii) studying the effects of the second Hunds 
rule in the orbital polarization theory pJ] , (iv) analyz- 
ing possible changes in the position of the Fermi level by 
changing the number of valence electrons pig , (v) us- 
ing the state tracking method JDjfl , and (iv) real space 
approach Q. 

In this paper we take a new view that the correla- 
tion effects within the d shell are important for the mag- 
netic anisotropy of 3d transition metals like Ni. These 
effects are not captured by the LDA but are described 
by Hubbard-like interactions presented in these systems 
and need to be treated by first principles methods JT|| . 

Another effect which has not been investigated in the 
context of magnetic anisotropy calculations is the non- 
collinear nature of intra-atomic magnetization |l?j| ]. It is 
expected to be important when spin-orbit coupling and 



correlation effects come into play together. In this arti- 
cle we show that when we include these new ingredients 
into the calculation we solve the long-standing problem 
of predicting the correct easy axis of Ni. 

We believe that the physics of transition metal com- 
pounds is intermediate between atomic limit where the 
localized d electrons are treated in the real space and fully 
itinerant limit when the electrons are described by band 
theory in k space. A many-body method incorporating 
these two important limits is the dynamical mean-field 
theory (DMFT) 0. The DMFT approach has been ex- 
tensively used to study model Hamiltonian of correlated 
electron systems in the weak, strong and intermediate 
coupling regimes. It has been very successful in describ- 
ing the physics of realistic systems, like the transition 
metal oxides and, therefore, is expected to treat properly 
the materials with d or f electrons. 

Electron-electron correlation matrix t/ 7l727374 = 
(miiri3 \vc \ S SlS2 S S3S4 for d orbitals is the quantity 

which takes strong correlations into account. This matrix 
can be expressed via Slater integrals , i = 0, 2, 4, 6 in 
the standard manner. The inclusion of this interaction 
generates self-energy E 7l72 (ia> n , k) on top of the one- 
electron spectra. Within DMFT it is approximated by 
momentum independent self-energy £ 7l72 (iw ra ). 

A central quantity of the dynamical mean-field theory 
is the one-electron on-site Green function 



G 7l72 (iw„) = J2k ( iuj n + A*)0 7l72 (fc) - # 7l72 (fc) 



Vjc 



(1) 



where 7? 7l72 (fc) is the one-electron Hamiltonian stan- 
dardly treatable within the LDA. Since the latter already 
includes the electron-electron interactions in some aver- 
aged way, we subtract the double counting term Vdc [@ ■ 
The use of realistic localized orbital representation such 
as linear muffin-tin orbitals |l9| ] leads us to include over- 
lap matrix 11 - Y2 (k) into the calculation. 

The DMFT reduces the problem to solving effec- 
tive impurity model where the correlated d orbitals are 
treated as an impurity level hybridized with the bath of 
conduction electrons. The role of hybridization is played 
by the so-called bath Green function defined as follows: 
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7172 (W") — G 7l72 (iUJ n ) + S 7l72 

(2) 

Solving this impurity model gives access to the self- 
energy S 7l72 {iu) n ) for the correlated electrons. The one- 
electron Green function (Q) is now modified with new 
X 7172 (iu> n )i which generates a new bath Green function. 
Therefore, the whole problem requires self-consistency. 

In this paper we confine ourselves to zero temperature 
and make an additional assumption on solving the im- 
purity model using the Hartree-Fock approximation. In 
this approximation the self-energy reduces to 



^7172 



7172 7374 
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where 



■-7172 



is the average occupation matrix for the cor- 



related orbitals. The off-diagonal elements of the occu- 
pancy matrix are not zero when spin-orbit coupling is 
included PQj . The latter can be implemented following 
the prescription of Andersen [fl9[ or more recent one by 
Pederson |2jj ]. 

In the Hartree-Fock limit the self-energy is frequency 
independent and real. The self-consistency condition of 
DMFT can be expressed in terms of the average occu- 
pation matrix: Having started from some n 7l72 we find 
S 7l72 according to (||). Fortunately, the computation of 
the on-site Green function ([!]) needs not to be performed. 
Since the self-energy is real, the new occupancies can 
be calculated from the eigenvectors of the one-electron 
Hamiltonians with S 7l72 — Vd c added to its dd block. The 
latter can be viewed as an orbital-dependent potential 
which has been introduced by the LDA+U method |lS{ ]. 

The LDA+U method has been very successful com- 
pared with experiments at zero temperature in ordered 
compounds. By establishing its equivalence to the static 
limit of the DMFT we see clearly that dynamical mean- 
field theory is a way of improving upon it, which is crucial 
for finite temperature properties. 

In this work we study the effect of the Slater parame- 
ter Fq that is the Hubbard on site interaction U on the 
magnetic anisotropy energy. The latter is calculated by 
taking the difference of two total energies with differ- 
ent directions of magnetization (MAE=E(111)-.E(001)). 
The total energies are obtained via fully self consistent 
solutions. Since the total energy calculation requires high 
precision, full potential LMTO method [j22j has been em- 
ployed. For the k space integration, we follow the analy- 
sis given by Trygg and co-workers Jll[ and use the spe- 
cial point method |2^] with a Gaussian broadening ]2^] 
of 15 mRy. The validity of this procedure has been 
tested in their work |yj. For convergence of the total 
energies within desired accuracy, about 15000 fc-points 
are needed. We used 28000 fc-points to reduce possible 
numerical noise. Our calculations include non-spherical 
terms of the charge density and potential both within 
the atomic spheres and in the interstitial region j2^] . All 
low-lying semi-core states are treated together with the 
valence states in a common Hamiltonian matrix in order 



to avoid unnecessary uncertainties. These calculations 
are spin polarized and assume the existence of long-range 
magnetic order. Spin-orbit coupling is implemented ac- 
cording to the suggestions by Andersen We also 
treat magnetization as a general vector field, which re- 
alizes non-collinear intra-atomic nature of this quantity. 
Such general magnetization scheme has been recently dis- 
cussed in ]l6| . 

We now discuss our calculated MAE. We first test our 
method in case of LDA (Z7 = 0). To compare with previ- 
ous calculations, we turn off the non-collinearity of mag- 
netization which makes it collinear with the quantization 
axis. The calculation gives correct orders of magnitude 
for both fee Ni and bec Fe but with the wrong easy axis 
for Ni, which is the same result as the previous result jll| . 
Turning on the non-collinearity results in a a larger value 
of the absolute value of the MAE (2.9 ^eV) for Ni but 
the easy axis predicted to be (001) which is still wrong. 
The magnitude of the experimental MAE of Ni is 2.8 neV 
aligned along (111) direction pH]. 

We now describe the effect of strong correlations, which 
is crucial in predicting the correct axis of Ni (see Fig. g). 
As U increases, the MAE of Ni smoothly increases until 
U reaches 2.5 eV and then smoothly decreases up to the 
value 3.8 fj,eV. Around U = 3.9 eV, the MAE decreases 
abruptly to negative value. Around U = 4.0 eV, the 
experimental order of magnitude and the correct easy 
axis (111) are restored. The change from the wrong easy 
axis to the correct easy axis occurs over the range of 
SU ~ 0.2eV, which is the order of spin-orbit coupling 
constant (~ O.leV). 

For Fe, the MAE decreases on increasing U to negative 
values, where the magnetization takes the wrong axis. 
From U = 2.7 eV, it increases back to the correct direc- 
tion of easy axis (positive MAE). Around U = 3.5 eV, it 
restores the correct easy axis and the experimental value 
of MAE is reproduced. 

It is remarkable that the values of U necessary to re- 
produce the correct magnetic anisotropy energy are very 
close to the values which are needed to describe photoe- 
mission spectra of these materials |2(|. This shows an 
internal consistency of our approach and emphasizes the 
importance of correlations. 

We find direct correlation between the dependency of 
the MAE as a function of U and the difference of mag- 
netic moments (AM = -(M(lll) - M(001)) behaving 
similarly. For Ni, the difference of magnetic moments is 
nearly U independent up to U = 3 eV. For large U, it 
smoothly decreases from the positive value to the nega- 
tive one. It also decreases rapidly around U = 3.9 eV 
in accord with MAE. For Fe, the difference is positive 
at U = Q. It decreases slightly to the negative values 
and then increases to the positive value over the range 
of U < 2.7 eV where MAE decreases. For larger U's, 
where MAE is coming back to positive value, its slope is 
significantly larger than that at smaller U's. 

This concurrent change of MAE and the difference of 
magnetic moments suggests why some previous attempts 
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based on force theorem |12|] failed in predicting the cor- 
rect easy axes. Force theorem replaces the difference of 
the total energies by the difference of one-electron ener- 
gies. In this approach, the contribution from the slight 
difference in magnetic moments does not appear and, 
therefore, is not counted in properly. Unfortunately, we 
could not find any experimental data of magnetic mo- 
ments to the desired precision (1CP 4 /xb) to compare with. 
We also have problems in reaching the convergence of the 
total energy with desired accuracy for large values of U 
in both Fe and Ni. 

We now present implications of our results on the cal- 
culated electronic structure for the case of Ni. One im- 
portant feature which emerges from the calculation is 
the absence of the X2 pocket (see Fig. ||). This has 
been predicted by LDA but has not been found exper- 
imentally ^7j. The band corresponding to the pocket 
is pushed down well below the Fermi level. This is ex- 
pected since correlation effects are more important for 
slower electrons and the velocity near the pocket is rather 
small. It turns out that the whole band is submerged un- 
der the Fermi level. 

There has been some suspicions that the incorrect po- 
sition of the X2 band within LDA was responsible for the 
incorrect prediction of the easy axis within this theory. 
Daalderop and coworkers removed the X2 pocket by 
increasing the number of valence electrons and found the 
correct easy axis. We therefore conclude that the absence 
of the pocket is one of the central elements in determining 
the magnetic anisotropy, and there is no need for any ad- 
hoc adjustment within a theory which takes into account 
the correlations. 

We now describe the effects originated from (near) de- 
generate states close to the Fermi surface. These have 
been of primary interest in past analytic studies p8| , p9[ . 
We will call such states degenerate Fermi surface crossing 
(DFSC) states. The contribution to MAE by non-DFSC 
states comes from the fourth order perturbation. Hence 
it is of the order of A . The energy splitting between 
DFSC states due to spin-orbit coupling is of the order 
of A because the contribution comes from the first order 
perturbation. Using linear approximation of the disper- 
sion relation e(fcA), the relevant volume in &;-space was 
found of the order A 3 . Thus, these DFSC states make 
contribution of the order of A 4 . However, there may be 
accidentally DFSC states appearing along a line on the 
Fermi surface, rather than at a point. We have found 
this case in our LDA calculation for Ni. Therefore the 
contribution of DFSC states is as important as the bulk 
non-DFSC states though the degeneracies occur only in 
small portion of the Brillouin zone. 

The importance of the DFSC states leads us to com- 
parative analysis of the LDA and LDA+U band struc- 
tures near the Fermi level. In LDA, five bands are cross- 
ing the Fermi level at nearly the same points along the 
TX direction. Two of the five bands are degenerate for 
the residual symmetry and the other three bands acciden- 



tally cross the Fermi surface at nearly the same points. 
There are two sp bands with spin up and spin down, 
respectively. The other three bands are dominated by 
d orbitals. In LDA+U, one of the d bands is pushed 
down below the Fermi surface. The other four bands 
are divided into two degenerate pieces at the Fermi level 
(see Fig. ^|): Two symmetry related degenerate bands 
and two near degenerate spi and spi bands. In LDA, 
we found that two bands are accidentally near degener- 
ate along the line on the Fermi surface within the plane 
TXL. One band is dominated by orbitals. The other 
is dominated by Sj orbitals near X and by orbitals 
off X. In LDA+U, these accidental DFSC states disap- 
pear (see Fig. ||). Instead, there is new two- fold DFSC 
states along TL direction, both of which are dominated 
by d[ orbitals. 

As we have seen, the on-site repulsion U reduces the 
number of DFSC states along TX direction while increas- 
ing that of DFSC states along TL direction. Based on the 
tight -binding model, Mori and coworkers |2!| have shown 
that DFSC states along TL direction result in the magne- 
tization aligned along (1, 1, 1) direction and DFSC states 
along TX direction result in the magnetization aligned 
along (0, 0, 1) direction. Since the strong correlation does 
precisely this, we conclude that disappearance of DFSC 
states along TX direction and their appearance along TL 
direction is another important element that determines 
the easy axis of Ni. 

Unlike LDA, we have found two extra very tiny L pock- 
ets in LDA+U (see Fig. ^|). Both of them are dominated 
by sp orbital with opposite spins. Being small, these ex- 
tra pockets may be artifacts of LDA+U. 

To conclude, we have demonstrated that it is possible 
to perform highly precise calculation of the total energy 
in order to obtain both the correct easy axes and the 
magnitudes of MAE for Fe and Ni. This has been accom- 
plished by including the strong correlation effect via the 
Hubbard-like on-site repulsion U and incorporating the 
non-collinear magnetization. In both Fe and Ni, the on- 
site U takes physically acceptable values consistent with 
the values known from atomic physics. The calculations 
performed are state of the art in what can currently be 
achieved for realistic treatments of correlated solids. Fur- 
ther studies should be devoted to improving the quality 
of the solution of the impurity model within DMFT and 
extending the calculation to finite temperatures. 
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FIG. 1. Experimental and calculated magnetocrystalline 
anisotropy energy MAE = £(111) - -5(001) (filled circle) and 
the difference of magnetic moment AM = Af(001) — M(lll) 
(square) for Fe (top) and Ni(bottom). Experimental MAEs 
are marked by arrows for Fe (1.4 [ieV) and Ni (—2.8 fieV). 




FIG. 2. Calculated Fermi Surface of Ni with the correla- 
tion effects taken into account. The solid and dotted lines 
correspond to majority and minority dominant spin carriers. 
Dominant orbital characters are expressed. Both experimen- 
tally confirmed X5 pocket and L neck can be seen. The X2 
pocket is missing, which is in agreement with experiments. 
There are two small L pockets, which has not been found by 
experiments. 
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